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CLT Abstract 
CD ' 

, ' This paper reports the result of our evaluation of the tenth-order QED correction to the lepton 

£Nj . g — 2 from Feynman diagrams which have sixth-order light-by- light-scattering subdiagrams, none 

> , 

of whose vertices couple to the external magnetic field. The gauge-invariant set of these diagrams, 

O 

called Set 11(e), consists of 180 vertex diagrams. In the case of the electron g — 2 (a e ), where 
the light-by-light subdiagram consists of the electron loop, the contribution to a e is found to be 

o , 

-1.344 9 (10) (a/vr) 5 . The contribution of the muon loop to a e is -0.000 465 (4) (a/vr) 5 . The 
contribution of the tau-lepton loop is about two orders of magnitudes smaller than that of the 



^ ■ muon loop and hence negligible. The sum of all of these contributions to a e is —1.345 (1) (a/V) 5 

h : 

We have also evaluated the contribution of Set 11(e) to the muon g — 2 (oA The contribution 
to from the electron loop is 3.265 (12) (a/V) 5 , while the contribution of the tau-lepton loop is 
—0.038 06 (13) (a/ir) 5 . The total contribution to a„, which is the sum of these two contributions 
and the mass-independent part of a e , is 1.882 (13) (a/ir) 5 . 

PACS numbers: 13.40.Em,14.60.Cd,14.60Ef,12.20.Ds 
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I. INTRODUCTION 



The anomalous magnetic moment g — 2 of the electron has provided the most stringent 
test of the validity of quantum electrodynamics, QED. The experimental value with the 
least uncertainty is the one obtained by the Harvard group in 2008 (a s fl 

a e (HV08) = 1 159 652 180.73 (28) x 1(T 12 . (1) 

To confront the prediction of the standard model with this measurement the hadronic con- 
tribution up to the order a 3 2|-|9|, the electroweak contribution up to the two-loop order 
ipj-12], and the QED radiative correction up to the eighth order must be taken into account 
13Hl5|. In order to match or exceed further improvement in the accuracy of the experimen- 
tal value, it is necessary to evaluate the tenth-order QED radiative correction to g — 2. To 
meet this challenge we launched several years ago the project to compute all 12672 Feynman 
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diagrams that contribute to the tenth-order a e 

The most difficult to evaluate is the gauge-invariant set, called Set V, which consists of 
6354 diagrams that have no virtual lepton loop. To deal with this set systematically we 



have developed an automatic code-generating algorithm GENCODEiV [16|, [18j. We now have 
FORTRAN codes for all diagrams of Set V generated by GENCODEiV. Numerical evaluation 
of these integrals is in progress at present. 

Meanwhile, we have also made steady progress in the evaluation of other types of tenth- 



order diagrams, and have published some of the results [17J, [19|, [20j . At the tenth order there 
appear five gauge-invariant sets of diagrams, called Set I(j), Set 11(e), Set 11(f), Set III(c) , 
and Set VI(j), which contain light-by- light-scattering subdiagram(s) internally, i.e., none of 



whose vertices is the external vertex 
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Uj. (See Fig.rU) 



Feynman diagrams containing a light-by-light-scattering subdiagram internally appear 
for the first time in the eighth-order QED correction to the lepton g — 2. Figure [2] shows 
the eighth-order self-energy diagrams with the fourth-order internal light-by-light-scattering 
subdiagrams. Vertex diagrams relevant to lepton g—2 can be obtained by inserting a single 
external QED vertex into one of open lepton lines labeled 1, 2, 3 of individual diagrams of 
Figure [3 

The diagrams of Set I(j) are those involving two fourth-order light-by-light-scattering 
subdiagrams, both internal, which have been evaluated and published recently [3]. The 
diagrams of other subsets are obtained by adding 0(a) correction to those of Fig. [2j The 
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I(j) 11(e) 11(f) III(c) Via) 



FIG. 1: Representative diagrams of gauge- invariant sets which contain an internal light-by-light- 
scattering subdiagram. 

Set 11(f) consists of diagrams obtained by inserting a second-order vacuum-polarization loop 
into one of internal photon lines of the diagrams of Fig. [2] in all possible ways. They have 
been evaluated by a simple modification of the FORTRAN codes developed previously for 
the eighth-order work. The result was published in Ref. The diagrams of Set III(c) 

are obtained by attaching a virtual photon line to the open lepton path of the individual 
diagrams of Fig. [2] in all possible ways. Evaluation of this set is in progress. The diagrams of 
Set VI (j) contain two light-by-light-scattering subdiagrams, one of which is internal, while 
the other is external. The numerical result of the Set VI (j) was published in Ref. [17] . 

The diagrams of Set 11(e) are obtained by attaching both ends of a virtual photon line 
to the lepton loop of the individual diagrams of Fig. [2] in all possible ways, forming the 
sixth-order internal light-by-light-scattering subdiagram. This paper reports the result of 
our work on Set 11(e), which consists of 180 Feynman vertex diagrams. 

The paper is organized as follows. Section [IT] describes the strategy we have adopted 
for the numerical study. The renormalization is set up so that ultraviolet divergences can 
be subtracted away without introducing spurious infrared divergence. Section II I II gives the 
results of our numerical work which covers the contributions of all diagrams of Set 11(e) to 
a e and a M . Section HYJ is devoted to the discussion and summary. 



II. COMPUTATIONAL PROCEDURE 



This section describes the strategy for computing the diagrams of Set 11(e). We denote 
the contribution of Set 11(e) to the magnetic moment of the lepton I induced by the virtual 
loop of lepton V as 

o,(n(e),0 = A(n(e),o(f) 8 . (2) 
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FIG. 2: Eighth-order self-energy Feynman diagrams LL8 that contain the fourth-order light-by- 
light-scattering subdiagram internally. In Ref. [2^], A, B and C are called "LLJ", "LLL" and 
"LLK" , respectively. The vertex diagrams are obtained by inserting a single QED vertex into one 
of the lepton lines 1, 2 or 3 of G (G = A, B, C). 
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where the lower case "a" includes the factor f-J while the upper case "A" does not. Recall 
that ai (11(e), /' = /) and Ai (11(e), I' = I) are independent of I and called mass- independent 
contributions. Ai (11(e), /' 7^ I) depends only on the mass ratio myjmi. We use the values 



found in Ref. 



2l| for lepton masses. 



As explained in Sec. HI the diagrams of Set 11(e) are obtained by attaching an internal 
photon line to the lepton loop of the individual diagrams of Fig. [2] in all possible ways. Let 
us denote the diagrams of Set 11(e) as Gij by specifying (i) the base eighth-order diagram G, 
where G is one of A, B, or C of Fig. |2l and (ii) a pair (i, j) of lepton lines of the closed loop 
(4 < i < j < 7) to which an additional internal photon line is attached. For instance, the 
insertion of two QED vertices into the middle of the lines 4 and 7 of the diagram A of Fig. [2] 
and the introduction of a virtual photon line which connects these vertices produces the 
tenth-order diagram called A47. Representative diagrams of Set 11(e) are shown in Figure 
The charge conjugation and time reversal symmetries of QED help us to reduce the number 
of independent amplitudes. For instance, A67 gives the same contribution to g — 2 as AA7. 
Recall also that the diagram in which the lepton loop runs in the opposite direction gives the 
same contribution as the original one. In this way, we obtain a complete set of independent 
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diagrams: 

AU [4], A55 [2], A77 [2], A47 [4], A45 [4], A46 [2], A57 [2], 
544 [4], £55 [4], B47 [4], 545 [4], 546 [2], 557 [2], 

C44 [4] , C55 [2] , C77 [2] , CA7 [4] , C45 [4] , C46 [2] , C57 [2] , (3) 

where the number in the brackets accounts for the symmetry factor for each diagram as well 
as the doubling due to two directions that a lepton loop takes. 

Thus far no one has succeeded in evaluating the diagrams of Set 11(e) a naly tically. We 
resort to the numerical means utilizing the parametric integral formulation [m, [rsl. 22 . 23]. 
The evaluation of g—2 can be simplified significantly by focusing on the quantity associated 
with the self-energy diagram Gij, such as the magnetic moment amplitude Moij, using 
the Ward-Takahashi identity which relates the regularized self-energy function J^Gijip) of 
the diagram Gij to the sum AaijiP, <?) of the contributions from the regularized vertex 
diagrams obtained by inserting a QED vertex into Gij in all possible ways [24]. 

The next step is to renormalize the integrals on the computer, which we carry out by 
subtractive renormalization. Since the bare amplitudes of individual diagrams have different 
structures of UV singularities, the numerical subtraction of UV singularities must be carried 
out for each diagram separately. Our aim is to construct subtraction terms that (i) share 
the same UV singularity as the integrand of the bare amplitude in the common Feynman 
parameter space, and (ii) do not introduce spurious IR singularities. 




1 2 3 1 2 3 1 2 3 

A47 AU A57 



FIG. 3: Representative diagrams of Set 11(e). ^147 and ^444 involve a second-order vertex subdia- 
gram and a second-order self-energy subdiagram, respectively. 
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The second point is not a trivial requirement. For instance, the usual on-shell second- 
order vertex renormalization constant contains an IR divergence. In general the subtraction 
term constructed under the usual on-shell renormalization condition introduces an IR sin- 
gularity that is not present in the bare amplitude. To avoid this problem we perform the 
renormalization in two steps. The first step is to construct the UV-finite amplitude AMaj in 
which only the UV-divergent part of the corresponding on-shell vertex (or self-energy) term 
is subtracted, leaving out the UV-finite piece unsubtracted. We call this step an intermedi- 
ate renormalization. The second step is to carry out the finite residual renormalization to 
account for the difference between the intermediate renormalization and the usual on-shell 
renormalization. The IR-divergent parts of the usual on-shell renormalization constants 
appear in the second step but cancel out when summed over the entire gauge-invariant set. 

The subtraction terms of AMoij are constructed as follows. The UV singularities as- 
sociated with the second-order vertex and self-energy subdiagrams are subtracted via K- 
operation, retaining the Feynman cut-off until UV divergences cancel out by renormaliza- 
tion [23[. The UV singularities of the light-by-light-scattering (/-/) loops are subtracted while 
maintaining the Pauli-Villars regularization in order to avoid dealing with divergent hence 
undefined quantities. The Pauli-Villars mass is sent to infinity only after renormalization is 
carried out. Note that not only the usual on-shell renormalization but also the intermediate 
renormalization are defined on the mass-shell insofar as it is IR-safe. To avoid confusion let 
us call the usual on-shell renormalization as the full renormalization henceforth. 

Let us illustrate our renormalization procedure taking A47 of Fig. [3] as an example. 
A47 has four UV-divergent subdiagrams which can be identified by the sets of lepton lines 
involved: 

S x = {1, 2, 4, 5, 6, 7, 8, 9} , S 2 = {2, 3, 4, 5, 6, 7, 8, 9} , 

S 3 = {4, 5, 6, 7, 8, 9} , S 4 = {8, 9} . (4) 

Both subdiagrams Si and 5*2 are the eighth-order vertex subdiagrams, S3 is the sixth-order 
/-/ subdiagram, and S4 is the second-order vertex subdiagram. Each UV subtraction term of 
AMcij is associated with a Zimmermann's forest. A47 has 11 normal forests. Let C$ denote 
the operator which extracts the full renormalization constant of the subset S from Maij, 
and let Kg denote the operator which extracts the UV singularity of the subset S by the 
intermediate renormalization defined by the K operation, respectively. Then the UV-finite 
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quantity AM A47 is denned by 

A.U ,,v = M A47 - C Sl M A47 - Cs,.W,.;v - C Sa M A47 - K Si M A47 
+K 54 C 53 M A4 7 + K 54 C 5l M A4 7 + Ks 4 C S2 M A47 
+C S3 C Sl M A47 + C S:i C S2 M A47 

-K S4 Cs 3 C Sl M A47 - K S4 C S3 C S2 M A47 . (5) 

Expression of AMaij for other Gij in Eq. ([3]) can be written down similarly. 

It is by definition that all subtraction terms on the right-hand side of Eq. (jSj) are fac- 
torizable. For instance, the operator C Sl acting on M^ 47 produces the product of L Sl and 
M 2 : 

C 5l M A47 = L Sl M 2l (6) 

where L$ 1 is the full eighth-order vertex renormalization constant of the diagram that con- 
tains the sixth-order light-by-light-scattering subdiagram, and M 2 = a 2 = | is the second- 
order lepton g— 2. Of course this equation is meaningless unless it is regularized. L,<?, and M 2 
can be expressed as regularized integrals in the parametric integral formulation [22[ on two 
separate Feynman parameter spaces with constraints J^i-aii lmeses x i = ^1 Sfc=3 dVk = 1> 
where yd is the Feynman parameter associated with the photon d. 



A manipulation similar to that of Sec. Ill D of Ref. [16( expresses Ls 1 M 2 as an integral 
over the single Feynman parameter space with YljeAu z i = ^- With this form of integrand 
of Ls 1 M 2 the pointwise subtraction of the overall UV divergence of M^ 47 residing in the sub- 
diagram S\ can be achieved. [Actually, Cs 1 Ma47 still has divergences from the subdiagrams 

53 and S 4 which must be subtracted by other terms of Eq. (J5]).] 

The K-operator, Ks 4 , acts on the regularized integrand J(z) of Ma 47 directly and produces 
a function Js 4 (z) that possesses the same UV singularity associated with the subdiagram 

5 4 as J{z). By definition K-operation also has the factorization property. For instance, the 
operator K$ 4 acting on the regularized M A47 produces the factorized result 

K S4 M A47 = L 2 v M 8A , (7) 



where L 2 Y is the UV-divergent part of the regularized second-order on-she 



malization constant L 2 and does not include the IR-divergent part of L 2 [23j. M 8A denotes 



1 vertex renor- 



the amplitude of the magnetic moment from the eighth-order diagram A of Fig. [2j The 



regularization mass must be sent to infinity after K$ 4 M AA7 is combined with M447. The 
difference of L 2 y and L 2 is accounted for at the stage of the residual renormalization. 
The subtraction term Cs 3 Ma47 can be written (somewhat) symbolically as 

Cs 3 Ma47 = n(0,0,0,0)M 4 , (8) 
Here II is a short-hand form of the sixth-order I- 1 sub diagram defined by 

^K\^u{k a , kb, k c , kd)\ka=k b =k c =k d =Oi (9) 

where k a , etc., are the momenta carried by the photon line a, etc., and M 4 is obtained 
from M AA7 by shrinking the /-/ loop of S3 to a point. The UV divergence of M AA7 arising 
from the subdiagram S 3 is cancelled by the term Cs 3 Ma47 of Eq. (jSJ), which results in full 
renormalization of the S3 divergence. Actually, M447 contains another UV divergence arising 
from S4 which we subtract by the operator Kg A . The complete removal of UV divergences 
arising from 5*3 and S4 is achieved by the combination 

M AA7 - C Ss M AA7 - K S4 M AA7 + K S4 C S3 M AA7 . (10) 

When the contributions of all diagrams of Set 11(e) listed in Eq. (jHJ) are put together, 
11(0, 0, 0, 0) from all diagrams cancel out and we obtain a simple result 

A(H(e), I') = AM ^J ~ 4 { ( L * - ^2 V ) + {B% ~ 5 2 UV ) } A(LL8, I') 

GiieEq.|3j 

= Yl AM Gi3 --4A5 2 x A(LL8, I'), 

GiieEq.|3j 

AB 2 = {(L 2 - LH + (B 2 - B^)} = 1. (11) 

Now, at last, we can send the regulator mass to infinity. B 2 is the full second-order wave 
function renormalization constant and B^ is the UV-divergent part of B 2 defined by the K 
operation. The IR divergence of [L 2 — L 2 y ) cancels that of (B 2 — B 2 V ) leaving a finite term 



as expected. (Note that B + L = while B% v + L^ v is finite but not zero. See Ref. [24] 
for the exact definitions of B 2 V and L% .) Eq. (iTTj) shows that the residual renormalization 
term is proportional to the eighth-order contribution v4/(LL8, I') to the anomalous magnetic 
moment of the lepton I from the diagrams of Fig. [2j in which loops are given by lepton I'. 



The numerical study in Ref. 25] has provided an accurate value for the mass-independent 
contribution v4 e (LL8, e) 

A e (LL8, e) = -0.990 72 (10) . (12) 
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In addition the paper [25J reports the electron-loop contribution to a M 



A M (LL8, e) 



-4.432 43 (58) . 



(13) 



We have also evaluated the muon loop contribution v4 e (LL8, fi) to the electron g — 2, and 
the tau-lepton loop contribution ^(LLS, r) to the muon g — 2 needed for this work 



v4 e (LL8, fj) 
A M (LL8, r) 



-0.000 177 8 (12), 
-0.015 868 (37). 



(14) 
(15) 



The remaining task is to evaluate every AMg^ in various combination of the external and 
internal leptons. 



III. NUMERICAL RESULTS OF AM Gij 

FORTRAN codes for AMcij are rather complicated and not easy to obtain. In order 
to facilitate this problem we adapted the automating code gencodeA specifically for the 
Set 11(e) which generates the integrands of AMaj as FORTRAN- formatted source programs. 



(See Refs. [16j, [18( for the details of automation.) Two independent sets of automating codes 
together with another set of manually-produced codes were constructed to confirm their 
validity. 

The integral for the diagram Gii, i.e, the one containing the second-order self-energy 
subdiagram, was found to exhibit worse convergence than the others. In order to alleviate 
this problem, we modify the integrand in the following way. For instance, in the diagram 
AA4 in Fig. [3j which contains a second-order self-energy subdiagram, the integrand of AM^ 
depends on the Feynman parameters Z4, z 8 only through the combination z 48 = (24 + z 8 ). 
Thus, the number of independent variables is reduced from 12 to 11. This seems to improve 
somewhat the convergence of iteration procedure. 

The integration of AM^ is carried out with the help of the adaptive-iterative Monte- 



Carlo integration routine VEGAS 26] on the massively parallel computer, RIKEN Integrated 
Cluster of Clusters (RICC). The number of sampling points for each iteration is 10 8 for all 
diagrams with the second-order self-energy subdiagrams and 2 x 10 s for all others. 
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TABLE I: Numerical results for mass-independent AMaj from diagrams Gij in Set 11(e). Full 
symmetry factors are included for the individual values. The number of sampling points for each 
iteration is 10 8 for all diagrams with the second-order self-energy subdiagrams and 2 x 10 8 for all 
others. The first, second, third and fourth columns show the name of the diagram, the value of 
integrand and its uncertainty, the x 2 value of VEGAS integration, and the number of iterations of 
VEGAS integration. If \ 2 is very close to 1, then the numerical integration by VEGAS is reliable. 



Gij 


AMcij (uncertainty) 


x 2 


iterations 


AAA 


5.397 41 (33) 


1.04 


2420 


A55 


2.796 88 (23) 


1.09 


1250 


A77 


2.422 84 (20) 


1.08 


1210 


A47 


0.100 26 (15) 


1.02 


1120 


A45 


-1.559 22 (16) 


0.98 


1150 


A46 


1.180 76 (12) 


1.10 


1140 


A57 


1.653 245 (93) 


1.09 


800 


BU 


-4.440 95 (34) 


1.06 


2660 


555 


-4.741 06 (33) 


0.99 


2500 


547 


1.725 96 (17) 


1.05 


990 


545 


2.521 96 (17) 


1.07 


1070 


546 


-0.349 57 (11) 


1.04 


1040 


557 


-2.254 206 (97) 


1.10 


790 


(744 


-5.054 64 (34) 


1.03 


2570 


C55 


-2.398 68 (24) 


1.09 


1670 


C77 


-2.431 20 (22) 


1.06 


1460 


C47 


1.574 62 (16) 


1.06 


990 


C45 


1.821 84 (17) 


1.00 


1200 


cm 


-1.881 49 (12) 


1.03 


1280 


C57 


-0.401 777 (91) 


1.03 


710 


sum 


-4.317 02 (94) 
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A. electron g— 2 



The results of numerical integration of AMaij for the mass-independent contribution to 
the lepton g—2 are presented in Table [H Following Eq. (fTTj) the last line of Table U together 
with the value (I12p for A e (LL8, e) yields the mass-independent contribution to g — 2 from 
Set 11(e) diagrams 

A, (11(e), I) = -1.344 86 (99), (16) 

where I = e,/z, or r. Recall that the actual contribution to g — 2, a e (II(e), e), is A e (11(e), e) 
times the factor (^) 5 . 

The electron g — 2 also receives the Set 11(e) contribution induced by the virtual muon 
loop. To see its numerical significance, the computation of AMaj for the muon loop is 
explicitly performed. The results are shown in Table [XTJ Putting together the last line of 
this table and the value ( TT4"1) of A e (LL8, /i) we obtain the muon loop contribution to the 
electron g — 2 from Set 11(e) diagrams 

A e (ll(e), y) = -0.000 465 (4). (17) 

The size of this contribution is less than the numerical uncertainty of the electron-loop 
contribution A e (11(e), e) given in Eq. ( 116ft . Since the tau- lepton loop contribution is expected 
to be about two-orders of magnitudes smaller than the muon loop contribution and hence 
negligible, we present the sum of Eqs. (JTBjl and (|T7|) as our current best value for the 
contribution to the electron g — 2 from Set 11(e) diagrams 

a e (II(e)) = -1.345 (1) ^ . (18) 
B. muon g — 2 

The main contribution of Set 11(e) to the muon g — 2 arises from the diagrams each of 
which is induced by an electron loop. We present the numerical result of AMaj f° r the 
electron- loop contribution in Table IIIIl This table shows that the sum of AMaj is an order 
of magnitude larger than that of the mass-independent AMaj m Table [H However, as is 
seen from Eq. ([TBI for ^(1X8, e), partial cancellation takes place between the first term 
and the second term [— 4AB 2 x A M (LL8, e)]. As a consequence, we have 

A„(II(e), e) = 3.265 (12) . (19) 
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TABLE II: Numerical results for AMcij of electron g— 2 from diagrams Gij in set 11(e) in each of 
which muon induces light-by-light scattering. Full symmetry factors are included in the individual 
values. The number of sampling points for each iteration is 10 8 for all diagrams with the second- 
order self-energy subdiagrams and 2 x 10 8 for all others. 



Gij AMQij (uncertainty) xlO 3 % 2 iterations 



,444 


11.424 89 (49) 


1.19 


1360 


,455 


5.822 60 (36) 


1.29 


640 


A77 


6.167 18 (33) 


1.23 


640 


A47 


3.674 56 (28) 


1.06 


320 


A45 


0.715 82 (275) 


1.14 


320 


A46 


3.296 85 (25) 


1.18 


400 


A57 


3.686 56 (20) 


1.55 


240 


544 


-6.068 57 (45) 


1.02 


1280 


555 


-7.947 942 (44) 


0.98 


1280 


547 


-1.252 29 (29) 


1.16 


240 


545 


0.775 54 (25) 


1.04 


320 


546 


-0.389 63 (26) 


0.91 


240 


557 


-3.590 13 (20) 


1.54 


240 


C44 


-7.361 42 (46) 


1.14 


1120 


C55 


-3.125 44 (31) 


1.17 


720 


C77 


-2.927 29 (34) 


1.33 


480 


C47 


0.213 13 (28) 


0.95 


240 


C45 


-0.306 03 (30) 


0.91 


240 


C46 


-3.423 30 (26) 


1.29 


400 


C57 


-0.383 70 (17) 


1.25 


240 


sum 


-0.998 6 (14) 
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TABLE III: Numerical results for AMoij of muon g — 2 from diagrams Gij in set 11(e) in each of 
which electron induces light-by-light scattering. Full symmetry factors included in the individual 



values. The number of sampling points for each iteration is 


1 n& fr\T* oil rliQf 
1U IOI d.11 Uld^ 


5 iams wiLii ine seconci- 


order self-enerj 


IY subdiagrams and 2 x 10 8 for all others. 








AMaj (uncertainty) 


x 2 


iterations 


AAA 


21.914 7 


(35) 


1.02 


2180 


4 55 


10.747 4 


(21) 


0.99 


1540 


477 


10.438 3 


(20) 


0.99 


1300 


447 


6.166 8 


(27) 


1.01 


1370 


44^ 


7.843 7 


(28) 


1.04 


1370 


AA(\ 


-13.679 5 


(14) 


0.98 


1050 


A ^7 


— 14.181 8 


(11) 


1.09 


810 


RAA 


-25.919 5 


(40) 


1.02 


3050 


uoo 


-25.634 7 


(37) 


1.03 


3050 


R47 


39.794 5 


(28) 


1.08 


1760 


R4^ 


41.011 


(28) 


1.11 


1810 


R4fi 


-16.936 4 


(15) 


0.99 


1130 


R^7 


-22.029 9 


(11) 


1.23 


840 


HA A 


-41.123 9 


(42) 


0.99 


3010 


C55 


-20.500 5 


(29) 


1.02 


1770 


C77 


-20.929 8 


(29) 


1.05 


1610 


C47 


46.365 7 


(30) 


1.05 


1810 


C45 


47.994 


(30) 


1.04 


1890 


C46 


-22.236 3 


(15) 


1.06 


1130 


C57 


-19.136 5 


(12) 


1.08 


810 



sum -10.032 (12) 
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TABLE IV: Numerical results for AMaj of muon g — 2 from diagrams Gij in set 11(e) in each 
of which tau-lepton induces light-by-light scattering. Full symmetry factors are included in the 
individual values. The number of sampling points for each iteration is 10 s for all diagrams with 
the second-order self-energy subdiagrams and 2 x 10 s for all others. 



Gij AMq,ij (uncertainty) \ 2 iterations 



,444 


0.422 177 (18) 


1.09 


1600 


,455 


0.215 057 (18) 


1.19 


480 


A77 


0.219 086 (19) 


1.18 


320 


A47 


0.112 235 (14) 


0.96 


240 


A45 


-0.008 963 (17) 


1.07 


160 


A46 


0.119 699 (16) 


1.37 


160 


A57 


0.138 883 (10) 


1.63 


160 


544 


-0.239 377 (18) 


1.02 


1400 


555 


-0.298 375 (18) 


0.99 


1280 


547 


-0.017 671 (16) 


1.11 


160 


545 


0.053 945 (16) 


1.16 


160 


546 


-0.014 195 (13) 


0.96 


160 


557 


-0.140 250 5 (80) 


1.43 


240 


C44 


-0.288 614 (18) 


1.04 


1280 


C55 


-0.125 456 (18) 


1.14 


400 


C77 


-0.120 731 (18) 


1.18 


320 


C47 


0.019 880 (15) 


1.13 


160 


C45 


0.011 773 (16) 


0.89 


160 


C46 


-0.129 852 (14) 


1.09 


240 


C57 


-0.014 913 9 (86) 


1.17 


160 


sum 


-0.085 66 (7) 
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Thus the electron- loop contribution (11(e), e) is not much larger than the muon loop 
contribution (11(e), fi) of Eq. f fT6|) . Since the sign of A^ (11(e), e) is opposite to that 
of A M (11(e), /i), we are curious about the role that the tau-lepton contribution A^ (11(e), r) 
might play. Table HVl shows the result of AM^ for the Set 11(e) contribution to the muon g— 2 
induced by the tau-lepton loop. Equation (TTTi) . together with the value (TTBT) of A M (LL8, r), 
gives 



which is two orders of magnitude smaller than v4 M (II(e), e) or v4 M (II(e), fi). Summing up 
Eqs. (TT^|) . (TTBT) and (12"UI) . the Set 11(e) contribution to the muon g — 2 is found to be 



IV. DISCUSSION AND SUMMARY 

In this paper, we computed the contribution to the lepton g — 2 from the tenth-order 
QED diagrams of Set 11(e) that contain the sixth-order light-by-light-scattering subdiagram 
internally. The use of Ward-Takahashi identity, as well as the symmetries of QED, reduces 
the computation of 180 Feynman diagrams to that of 20 integrals AMcij- The intermediate 
renormalization to define AMcij is chosen so that the UV divergence associated with the 
second-order self-energy or vertex subdiagram is subtracted away by K-operation. Meanwhile 
the UV divergence arising from the /-/ loop is subtracted by full renormalization. This leads 
to simplification of the final result as is seen in Eq. (TTTj) . 

The Set 11(e) contribution to the electron g — 2 is obtained by evaluating the electron 
and muon virtual effects. The result is given in Eq. (TT8"j) . The size is of the typical order of 
magnitude for the tenth-order. The numerical computation was carried out as accurately as 
possible with the available computer resources. 

The contribution to the muon g — 2 is obtained by evaluating the virtual effects of all 
leptons. The result is given in Eq. ([2"Tj) . The contribution of the electron loop to the muon 
g— 2 is not much larger than the muon loop contribution. 

We found that (11(e), e) in Eq. fTiTJl) involves partial cancellation between the sum of 
AMcij over all Gij in Set 11(e) and the residual renormalization term in Eq. fflTl) . In spite 



A,(n(e), r) 



0.038 06 (13) 



(20) 




(21) 
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of these problems we were able to obtain the result for a ^(1 1(e)) with the uncertainty less 
than 1% using the high performance computer system, RICC. 
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